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Description 

FIELD OF THE INVENTION 

5 [0001] The invention relates to electronic connnnunications, and more particularly to CDMA-based coding, transmis- 
sion, and decoding/synthesis methods and circuitry. 

BACKGROUND OF THE INVENTION 

10 [0002] Code division multiple access (COMA) coding has been extensively used in such applications as cellular and 
satellite communications. CDMA signals increase the spectrum required for the transmission at a particular data rate 
by modulating each data symbol with a spreading code having a rate larger than the data rate. The same spreading 
code Is used for each data symbol. Typically, the spreading code comprises of a few tens or a few hundreds of elements, 
called chips. To decrease the correlations among spreading codes assigned to different users, and thereby reduce the 

IS interference among different users, the data stream after spreading Is typically scrambled with a pseudo-noise (PN) 
code that Is generated serially and cyclically and has a larger period than the spreading code. Examples of such CDMA 
signal spreading are the schemes used by the IS-95/CDMA2000 and 3GPP systems. 

[0003] The 3GPP system employs a 3.84 MHz bandwidth (a chip duration of 260 ns), and the spreading code length 
applied to each data symbol may vary from 4 chips to 256 chips. The number of chips per symbol is the spreading 

20 factor (SF), and SF stays constant during a 10 ms duration frame. However, the spreading factor may change from 
frame to frame. Each frame is partitioned into 15 time slots (0.667 ms duration) with each time slot consisting of 2560 
chips, and thus the number of data symbols per slot ranges from 10 to 640 and per frame ranges from 150 to 9600. 
[0004] 3GPP allows fortransport channels with differing data rates andfonnats as specified by the Transport Format 
Control Indicator (TFCI) field, but 3GPP also allows for omission of the TFCI with instead only a (previously-negotiated) 

25 set of possible (candidate) transport block lengths known by the decoder. In this situation a decoder may determine 
the transport block length by decoding for all possible candidate block lengths and selecting one of the decodings. In 
particular, Figures 4-5 illustrate a 3GPP blind transport format detection (BTFD) method which uses the CRC bits 
(attached at the end of a transport block prior to convolutional encoding) as follows. First the received stream of symbols 
is Viterbi decoded up to the end of the shortest of the candidate block lengths. Viterbl decoding includes state metric 

30 updatings with add-compare-select (ACS) computations for each received symbol plus recording in the trace-back 
buffer which detected transition (0/1 input to the convolutional encoder) was used for each state metric update; this 
uses a buffer for current state metrics plus a buffer for the streams of inferred 0/1 inputs (the trace-back buffer). Next, 
the state metric for state 0 at the end of this shortest block is checked for reliability, and if reliable, then the trellis path 
ending at this state 0 is traced back to the beginning. The traced-back path provides a decoding for the shortest block 

35 length which can be error checked with its presumed CRC bits. Second, the Viterbi decoding is continued from the 
end of the shortest block to the end of the second shortest block. Again, the state metric for the state 0 at the end of 
the second shortest block Is checked for reliability, and if reliable, then the trellis path ending at this state 0 Is traced 
back to the beginning. This traced-back path provides a decoding for the second shortest block length which can be 
error checked with its presumed CRC bits. Similarly for the remainder of the block lengths; and the best (if more than 

40 one) of the error-free decodings is selected. This method uses a logarithm function in the reliability check of the state 
metric for the state 0. 

[0005] The known method has problems including large trace-back buffers and complex logarithm computation. 
[0006] Black et al, A 1-Gb/s, Four-State, Sliding Block Viterbi Decoder, 32 IEEE J. Solid State CIr. 797 (1997) de- 
scribes a Viterbi decoder using a sliding decoding block of length less than the entire block to decode. 

45 

SUMMARY OF THE INVENTION 

[0007] The present invention provides a low complexity implementation for blind forniat detection with a sliding win- 
dow Viterbi decoding. 

50 [0008] According to a first aspect of the invention there is provided a method of determining the length of a received 
encoded block, comprising: (a) receiving a sequence of encoded samples, said sequence including blocks of a length 
selected from a set of candidate lengths; (b) Viterbi decoding said sequence with a sliding window up to a first of said 
candidate lengths; (c) evaluating said decoding of step (b); (d) Viterbi decoding said sequence with a sliding window 
up to a next of said candidate lengths; (e) evaluating said decoding of step (d); (f) comparing the evaluating of step 

S5 (e) with prior evaluating; (g) repeating steps (d)-(f) to include all of said candidate lengths; and (h) determining a block 
length from said comparing of step (f). 

[0009] According to a second aspect of the invention there Is provided A method of decoding, comprising: (a) pro- 
viding a set of candidate block lengths; (b) providing a window length and a window slide length; (c) updating state 
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metrics with subtraction of minimum metric value of previous interval and starting at state 0 for a trellis using a sequence 
of received samples over an interval of length equal to said window length, said updating including storing correspond- 
ing transition bits; (d) tracing back from said updated state metrics with subtraction of minimum metric value of previous 
interval using said stored transition bits; (e) retaining bits of said tracing back from an end of said tracing back for a 

5 length equal to said window slide length; (f) updating said state metrics over a next interval of length equal to said 
window slide length, said updating including storing corresponding transition bits; (g) repeating steps (d)-(f) until said 
updating reaches a first of said set of candidate lengths; (h) tracing back from state 0 of said updated state metrics of 
step (g) using said stored transition bits, and combining with said retained bits for a first decoding up to said first of 
said set of candidate lengths, and evaluating said first decoding, and continuing updating said state metrics to the end 

10 of said next interval of length equal to said window length; (i) repeating steps (d)-(f) said updating reaches a next of 
said set of candidate lengths; (j) tracing back from state 0 of said updated state metrics of step (t) using said stored 
transition bits, and combining with said retained bits for a next decoding up to said next of said set of candidate lengths, 
and evaluating said next decoding, and continuing updating said state metrics to the end of said next interval of length 
equal to said window length; (k) repeating steps (i)-(j) until all of said set of candidate lengths has been traced back; 

IS and (I) comparing said evaluatings to select one of said set of candidate lengths. 

[0010] According to a third aspect of the invention there is provided a logarithm approximation method, comprising: 
(a) providing N-bit numbers c and d where N a positive integer and d is less than or equal to c; (b) extracting M-bit 
number m from said c where M is a positive Integer less than N, said m aligned to the most significant non-zero bit of 
c; (c) extracting M-bit number n from said d, said n aligned to said m; (d) piecewise lineariy approximating a logarithm 

20 of said m; (e) adjusting the approximation of step (d) by a table entry; (f) piecewise linearly approximating a logarithm 
of said n; (g)adjusting the approximation of step (f) by a table entry; and (h) approximating a logarithm of the quotient 
d/c by the difference of results of step (g) and step (e). 

[001 1 ] According to a fourth aspect of the invention there is provided a method of logarithm approximation , compris- 
ing: (a) providing N-bit positive integer c = x^X2 . . . x^^ with N a positive integer and x^ the most significant bit and x^^ 

25 the least significant bit; (b) providing N-bit positive integer d = yiya - yN with y^ the most significant bit and y^ the least 
significant bit and with d less than or equal to c; (c) extracting M-bit number m from said c, said m with most significant 
bit aligned to Xj of said c where x^, X2, X|(.^ are all 0 and is 1 where k is a positive integer with 1 ^ k < N; (d) 
extracting M-bit number n from said d, said n with most significant bit aligned to with said k from step (c); (e) when 
said m is in the range 2' to 2"+"* -1 for a positive integer I, define a number t = m2^"' + (M-i)2'^*i + table [m] where table 

30 [ ] is a lookup table; (f) when said n is in the range 2) to 2}^^ -1 for a positive integer], define a number u = n2M-i-J + (j- 
1)2 ^''^ + table[n]; and (g) taking t-u as an approximation of log(d/c). 
[0012] The preferred aspects provide more efficient blind transport fonnat detection. 

BRIEF DESCRIPTION OF THE DRAWINGS 

35 

[0013] The drawings are heuristic for clarity. 



45 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

1 . Overview 

[0014] Some preferred embodiment methods implement blind transport format detection (BTFD) as in 3GPP with a 
50 piecewise linear approximation of a logarithm for the reliability measure and other preferred embodiment methods use 
a windowed decoding. Figure 1a schematically shows a preferred embodiment windowing decoder with a reliability 
block which may include a preferred embodiment logarithm approximation as illustrated in Figure lb. 
[0015] Preferred embodiment systems Include cellular CDMA communication systems where several mobile users 
communicate with a base station, and the preferred embodiment BTFD methods and decoders are used. 
55 [0016] The preferred embodiment methods may be Implemented on programmable processors (such as DSP's) or 
in specialized circuits (ASIC's) or combinations such as system-on-a-chip with attendant analog-digital signal compo- 
nents. The methods could be stored programs in ROM with registers and on-board RAM for temporary data storage. 



40 



Figures 1a-1b are block diagrams of a preferred embodiment. 

Figure 2 illustrates encoders. 
Figure 3 shows a trellis butterfly. 

Figures 4-5 illustrate a blind transport fonnat detection method. 
Figures 6-9b show trellis paths. 
Figure 10 is a memory map. 
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2. Blind transport format detection 

[001 7] The first preferred embodiment blind transport format detection (BTFD) methods follow the methods of 3GPP 
BTFD using CRC but with a special approximation for the logarithm function; thus preliminarily consider 3GPP BTFD 
5 using CRC. 

[0018] 3GPP has two convolution coders with coding rates 1/2 and 1/3 with as illustrated in Figure 2 for shift register 
implementations. Each time increment provides an Input bit plus shifts the stored bits to generate the output bits (sym- 
bol); either Yq.Y^ for rate 1/2 or Yq, Y^, Y2 for rate 1/3. The pattern of the eight stored bits (the outputs of the eight 
delay units) defines the state of the shift register, and these 256 possible states are numbered 0 to 255 simply by 

10 Interpreting the bit pattern as a binary number with the LSB being the rightmost (and most recent) bit. Thus at each 
time interval, the MSB Is lost, the 7 remaining bits are shifted (multiplication by 2 of the 7-bit number), and a new LSB 
Input. That Is, the states 0x3X5X4X3X2X^X0 and 1x5X5X4X3X2X^X0 both transition to state XgX5X4X3X2X^XoO in the next time 
Interval with an input 0 and both transition to state X6X5X4X3X2X1X0I with an input 1 . More succinctly, if the 7-bit number 
X6X5X4X3X2X-1X0 is denoted n, then states n and n+128 both yield state 2n with an input 0 and both yield state 2n+1 with 

15 an input 1 . Graphically, the two start states, the two end states, and the four connecting transitions between them form 
a "butterfly"; see Figure 3. 

[0019] For encoding a block of bits, the initial state of a shift register of Figure 2 is eight Os (the state 0). And after 
encoding the block of bits, input a sequence of eight Os (tail bits) to force the shift register into state 0 for termination. 
[0020] All sequences of the possible 256 states and connecting transitions (0/1 inputs and corresponding Yo, Y^, 
(Y2) outputs) form the code trellis, and a particular sequence of states with connecting transitions fomns a path through 
the trellis. Thus a sequence of information bits input to the encoder (shift register) fomris a path through the trellis, and 
Viterbi decoding the received sequence of sets of soft symbols which correspond to modulation by Yq, Y^ for rate 1/2 
or Yq, Y^, Y2 for rate 1/3. Viterbi decoding estimates this path as the maximum likelihood path by computing a state 
metric for each of the 256 states and updating all of these state metrics upon each receipt of a symbol. More particularly, 
25 for the rate 1/3 encoder, let metric(i,n) be the state metric for state n at time i, then recalling states n and n+1 28 each 
have transitions to states 2n and 2n+1 . and presuming an AWGN transmission channel: 
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metric(i+1 , 2n) = 

max{metric(i,n)+2o[2Yo(n.2n)-1]+2i[2Yi(n,2n)-1]+22[2Y2(n,2n)-1], 
metric(i,n-i-128)+Zo [2Yo {n-i-128,2n)-13 +Zi[2Y^ (n+128,2n)- 
1]+Z2[2Y2(n+128,2n)-1]) 

metric (1+1 , 2n+1 ) = 
max{metric(l,n) + Zo[2Yo(n,2n+1)-1] + Zi[2Y^(n, 2n+1)-1] + 
Z2[2Y2(n,2n+1)-1], 
metric(i, n+128) + Zo[2Yo(n+128, 2n+1)-1] + ZJ2Yi (n+128,2n+1)- 
1] + Z2[2Y2(n+128, 2n+1)-1] } 



where Y|^(n, 2n) for k=0,1 ,2 are the three encoder output bits for a 0 input bit transition from state n to state 2n, Y|^ 
(n+128,2n) are the three output bits of the encoder for a 0 input transition from state n+128 to state 2n, and so forth. 
50 Note that [2Y^^- 1] converts the 0,1 of to -1,+1 of the modulation, and the {Zq, Z^, Z2} are the received detected 
(digitized) analog symbol corresponding to the transmitted output {Yo, Y^, Y2}. 

[0021 ] At each transition (time increment) for each state, store In a trace-back buffer one bit of Information to indicate 
which of the two entering transitions (prior states) was selected by the max{ }. In particular, at time i+1 for state 2n 
store a '0' when max{ ) came from the temn metric(i,n)+ Zo[2Yo(n,2n)-1]+ Zi[2Yi(n,2n)-1]+ Z2[2Y2(n,2n)-1] which cor- 
55 responds to an encoded bit 0, and store a *^' when the max came from the term metric(i,n+128)+ Zo[2Yo(n+128,2n)- 
1]+ Zi[2Yi(n+128,2n)-1]+ Z2[2Y2(n+128,2n)-1] which corresponds to an encoded bit 1 . This requires 256*N bits for a 
sequence of length N (N received {Zq, Z^, Z2)). At the end of the sequence (path) a trace is carried out from the final 
state 0 (forced by the termination with tail bits) backwards through the trellis to the initial state 0 as illustrated in Figure 
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6; this path (and the corresponding bits) is the Viterbi decoding. If all path selection information (a bit for each state at 
each transition) for a sequence of length 51 2 Is stored in the trace-back buffer, then the necessary buffer size is about 
128 kbits. 

[0022] As an alternative to explicitly transmitting TFCI. 3GPP allows for transmission once of a set of possible (can- 
5 didate) transport block lengths and relies on the receiver to use a blind transport format detection (BTFD) scheme to 
determine the actual block length. Figure 4 illustrates four possible block lengths (possible ending bit numbers n^j^^) 
and a block ending at ng^jj = 3. Figure 5 shows the flow of the BTFD which uses trellis trace back for each possible 
(candidate) block length together with the last bits as CRC check bits to check for errors for this block length. In order 
to limit the probability of the CRC missing errors, the BTFD method uses the reliability figure r{n^f^) in order to select 
10 a block length from candidate lengths. 

r ("end) = -10 logio ( [ao (riend) -amln ("end) 1 / [ ^niax ("end) '^mln (nendl ) 



survivor paths at the nend ending of the block, and (nend) 's the path-metric value of the 0 state at this block end. 
Note that by using the logarithm of the ratio, the reliability function is independent of the number of bits of accuracy of 
the state metrics. The method selects the correct, if any, block end n^^d using the following rule: 

20 1) Initialize the minimum reliability figure, r^,,,, to an available maximum value D. 

2) Conduct Viterbi decoding for one of the candidate block endings, n^n^j. 

3) Calculate the reliability figure r (n^nd) fo>' this candidate length. 

4) If r(ngnd) is greater than threshold D, return to step 2) and next candidate. 

5) Trace back from the 0 state at hgn^. 

25 6) Calculate CRC of the thus decoded bit sequence. 

7) If CRC Indicates error, return to step 2) for the next candidate. 

8) if r (n^nd) ^ ^miw return to 2) for the next candidate. 

9) Retain the sequence, update r^nin to r (n^,^), and return to step 2) for the next candidate. 

10) If there is no sequence retained by step 9), declare an error. 

30 11 ) Use the retained path from step 9) as the decoded sequence. Figure 5 is the flow diagram for this rule. 

[0023] Because the reliability figure r(ngp,jj) uses a logarithm, a look up table would be required to avoid introducing 
software or hardware to evaluate the logarithm function, and the lookup table would have size 8kbits having 1 k entries 
(10-bit data) with each entry 8 bits. This size lookup table is a problem. 
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3. First preferred embodiment blind transport format detections 

[0024] The first preferred embodiment blind transport fomnat detection methods reduce the complexity of a large 
lookup table or software implementation for implementation of the logarithm function used in the reliability function 

r (Hend) = »OgiO ( [^0 (^end) •^m\n ("end) ] ^ [^^^ (n^^^) -a^,^ (n^^^)] ) 
= -10 log^O ("end) -amin ("end) ) + 10 'OQlO (^max ("end) -amln ("end) ) 

by replacing the large lookup table with a 32-bit table plus a small amount of logic. Indeed, approximate each of the 
logarithms separately, noting that a constant added to both will cancel out. In particular, let 



55 so r(nend) = "10 log^o(^c) - 10 log-ioC - 10 logiod- Suppose a^^a^ and a^i^ are 10-bit numbers and let m be a 5-bit 
positive integer extracted from c such that the first non-zero bit of c corresponds to the MSB of m and the other 4 bits 
of m are the next 4 bits of c or Os if the c bits run out. Similarly, extract the 5-bit integer n from d in the same bit locations 
as m. Hence, to 5 significant bits, d/c = n/m and thus 
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^ ("end) = 10 'oQiqC - 10 logiod « 10 log^om - 10 log^o"- 

Note that the number of bits in a^ax> ^min> ^ ^® varied without affecting the 5-bit approxinnation; in particular, 
16-bit a^ax> ^miw ^ same definitions. Also, note that by definition m has an MSB equal to 1 (so 31 ^ 

m > 1 6) and m > n > 0. Also, if c=0 (so there is no non-zero bit), then a^^x = ^min = there is no reliability problem. 

[0025] Further, changing the base for the logarithm function only introduces a constant multiple but can make the 
lookup table simple. In particular, express the reliability with logarithms to base b: 



''("end) = 10 'OQio"^ - 10 log^on 
= logiob [10 logfam - 10 logfan] 

15 and the log^ob term can be absorbed into r{n^j^^) and D. Thus use a redefined reliability: 

r' ("end) = 10 logfam- 10 logbn 

20 A first preferred embodiment method picks b so that 10 log^ =16; that Is, take b = 2i^i-6 » 1 .54..., and so log^gb = 
0.1 875... Then the preferred embodiment approximates the first term 1 0 logjam of the reliability r* (n^^) by the variable 
t as follows 

25 if (m == 0) t = — 128; //this corresponds to 

else if (m < 2) t = table [m] ; 

30 else if (m < 4) t = table [m] + 8m; 



35 else if (m < 8) t = table [m] + 4m + 16; 

else if (ra < 16) t = table [m] + 2m + 32; 
else t = table [m] + m + 48; 

40 

where the array table[ ] has 32 entries which are either 0 or 1 as: 

table[32] = {0,0,0,1,0,1,1.1, 0,1,1.1,1,1,1,1, 
0.0.1,1,1.1,1,1, 1,1,1,1,1,1,1,0} 

Note that table[0] is not used (t = -128) and table[1] == 0 just reflects log1 = 0. 
[0026] Similarly, approximate the second term 10 log^^n by the variable u as follows: 
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if (n == 


0) 




u = -128; 




else if 


(n 


< 2) 


u - table [n] ; 




else if 


(n 


< 4) 


u = table [n] + 8n; 




else if 


(n 


< 8) 


u = table [n] + 4n + 


16; 


else if 


(n 


< 16) 


u = table [n] + 2n + 


32; 


else 






u = table [n] + n + 48; 





Consequently, approximate the redefined reliability r' (n^nd) 
15 [0027] For example, if a^ax = 1011010011. amm = 1001101001, and ao = 1010010010, then a^^^x - ^min = c = 
0001 101010 (=1 06), m = 1 1 01 0 (= 26), - a^^ = d = 00001 01 001 (=41 ), and n = 01 01 0 (= 1 0). Hence, t = table[26] 
+ 26 + 48 = 1 + 26 + 48 = 75 and u = table[1 0] + 2*1 0 + 32 = 1 + 20 + 32 = 53. Thus t-u = 75 - 53 = 22; so the 
approximated redefined reliability r* (n^^^) = 22. For comparison, the exact result would be 

20 

r(nend)= '•Ologioc- lOlogiod 
= 10logio106-10log^o41 
25 =20.253-16.128 

= 4.125 

and this agrees with the preferred embodiment approximation: 

30 

r(nend) = 'ogiob r' (n^nd) 
= 0.1875 r»(n^„d) 

= (0.1875)(22) 
= 4.125. 

40 [0028] The foregoing approximation is a piecewise linear approximation to the logarithm with the table as a correction 
term for the convexity of the logarithm. And Figure 1 a Illustrates (max, min ... block and reliability block) the computa- 
tions as first finding a^^jj, a^i^, and Qq; next, computing c and d (a subtraction circuit or programmed processor); next, 
finding m and n (a bit detection and extraction circuit or programmed processor); and then computing r' (n^^^) using 
multiplication, addition, and subtraction together with the lookup table. Figure lb shows the approximation functional 

45 blocks. 

4. Second preferred embodiments 

[0029] Second preferred embodiment methods are analogous to the first preferred embodiment methods but use 
50 other than 5-bit approximations in the logarithm. For example, with 4-bit approximation (so the lookup table has 16 
bits), the integers m and n are 4 bits and the MSB of m is again the first nonzero bit of c. and n is the 4 bits at the same 
locations in d. In this case the logarithm approximation would approximate lOlogm by the variable t as follows: 
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if (m == 


0) 






t 


= -128; 






else if 


(m 


< 


2) 


t 


= table [m] ; 






else if 


(m 


< 


4) 


t 


= table [m] + 


8m; 




else if 


(m 


< 


8) 


t 


= table [m] + 


4m + 


16; 


else 








t 


= table [m] + 


2m + 


32; 



where the array table[ ] has 1 6 entries which are Just the first haif of the first preferred embodiment tabie: table[16] = 
{0,0.0.1,0.1.1.1, 0,1,1.1.1,1.1,1} 

[0030] For the same exannple as in the first preferred embodiment, c (=1 06) and d (=41) remain the same, but the 
15 4-bit m = 13 and 4-bit n = 5. Thus t = table[13] + 2*13 + 32 = 59, u = table[5] + 4*5 + 16 = 37, and t-u = 22, the same 
as with the first preferred embodiment; so the approximation is the same. 

[0031] Further, the 4-bit preferred embodiment could also simplify the lookup table by changing the base for the 
logarithm from b = 1 .54... to b = (1 .54.. .)2 which will make 10iogt,2 = 8 and change the approximation to: 



20 





if (m == 


0) 






t = 


-128; 






else if 


(m 


< 


2) 


t = 


table [m] ; 




25 


else if 


(m 


< 


4) 


t = 


table [m] + 


4m; 




else if 


(m 


< 


8) 


t = 


table [m] + 


2m + 8; 




else 








t = 


table [m] + 


m + 16; 



30 

and similarly for u In terms of n. Because this essentially just divided the prior by 2, the lookup table changes to table 
[1 6] = {0,0,0,1 ,0,1 .1 ,0.0.0,0,1 .1.1 ,0,0} 

[0032] With this new b, the redefined reliability function r" {riend) = leads to: 

r ("end) = "Ogiob ^ ("end) 

= 0.375 r 

40 And again using the same example, c (=106), d (=41), 4-bit m = 13. and 4-bit n = 5 all remain the same, but t = table 
[13] + 13 + 16 = 30. u = table[5) + 2*5 + 8 = 19, and t-u = 11; so r (n^nd) = 0.375 r" (n^^^) = (0.375) (11) = 4.125. the 
same as before. In short, the approximation and table can be adjusted to the number of bits of resolution desired. 
[0033] Analogously, the number of bits in m and n can be greater than 5; this stili allows simple logic with an expansion 
of the lookup table and changes in the base for the logarithm. For example, positive m and n with M bits (in the range 

45 0 to 2^-1) and the logarithm base such that 10logb2 = 2^-1, if n is in the range 2i to 2J+'' -1 for a positive integer], take 
u = n2^-''-i + (j-1 )2'^-i 4- tabie[n] as the approximation of 1 0logbn where tab!e[ ] is a table of size 2^ (although the special 
cases of table[01 is irrelevant because u is taken as large negative and table[1] = 0 due to log^l = 0 for every base b) . 
Of course, for greater resolution the lookup table could be expanded to 2-bit or larger entries. 

so 5. Trace-back windowing preferred embodiments 

[0034] The trace-back windowing preferred embodiments reduce the complexity of foregoing preferred embodiment 
BTFD methods which use a Viterbi decoding for each of the candidate block sizes to select the correction one. In 
particular, the trace-back windowing preferred embodiments use sliding windows for the trace backs; this reduces the 
S5 trace-back buffer size. Indeed, consider a sliding window of size W. If W is sufficiently large, then a path traced back 
from the state on the end edge of the sliding window with the largest state metric will converge to the path obtained 
by the normal trace back from state 0 at the block end. Figure 7 illustrates this with trace-back paths from window end 
edges at W-1 , W+d-1 , and W+2d-1 plus the nomriai trace-back path. 
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[0035] The preferred embodiment sliding window decoding method with window length W and window slides of length 
dis: 

1) Start state metric calculation from state 0 at time 0; this is the usual. 
5 2) Continue state metric calculations up to the first window end edge at time W-1 . 

3) Carry out trace-back from the state with the largest state metric at time W-1 back to time 0. 

4) Retain the last d bits (corresponding to times 0, 1, d-1) out of the W decoded bits of step 3) and discard the 
remaining W-d bits. 

5) Continue the state metric calculations up to the next window end edge at time W+d-1 . 

10 6) Carry out trace-back from the state with the largest state metric at time W+d-1 back to time d (i.e., the front 

edge of this next window). 

7) Retain the last d bits (corresponding to times d, d+1 , 2d-1 ) out of the W decoded bits from step 6) and store 
them with the prior d bits from step 4); discard the remaining W-d bits. 

8) Repeat steps 5) - 7) until the state metric calculations reach the end of the input block. 

15 9) Carry out trace-back from the state 0 at the block end back to the front edge of the last window (which included 

the block end). 

10) Retain all decoded bits from step 9) and store together with prior retained and stored bits from step 4) and all 
repeats of step 7). 

20 Figure 7 illustrates the windows and trace-backs with k denoting the number of window slides to cover the entire block. 
[0036] The decoding quality (i.e., comparable to normal trace-back) of the preferred embodiment sliding window 
method depends upon the window size W and the sliding distance d (or window overlap size W-d). For example, Figure 
8 shows the case of a block of size 120, a window of size 80, and a slide distance of 16 (overlap size 64); so four 
windows cover the block. In general, a trace-back of length at least the survivor path length L = 5 log2M, where IVI is 

25 the number of states, converges to the correct path with very high probability. Thus a trace back of at least 40 time 
intervals (received symbols) should suffice. 

[0037] 1 6 bits are obtained in each of first three trace-backs of 80-bit window sequences and are illustrated as length 
1 6 paths Pq, Pi, and Pg in Figure 8. The final 72-bit trace back P3 starts from state 0 at the block end and traces back 
to the front of window W3. The concatenation of paths Pq. P2. and P3 provides the whole decoded sequence. 

30 [0038] In the sliding window method, once a trace-back is completed for a window of size W, the d decoded (retained) 
bits will not be influenced by the subsequent decoding process. Decoding of these d bits is complete and fixed. There- 
fore, the decoded bits for different code lengths vary only within the decoded bits from the last window whose trace 
back starts at state 0 of the end of the block. In the example of Figure 8, paths Pq, Pi, and Pg are common to all 
decoded bits with block length 120 plus blocks larger than 120. But P3 is valid only for block length of 120. Presuming 

35 the block length 120 example of Figure 8 is one of 16 candidate block lengths (denoted Cq, C^, .... C^g) and that 16 
bits can be retrieved by a single trace back window, then further windows would be used. 

[0039] To decode the input sequence repeatedly, each decoding with one of the 1 6 different block lengths, then the 
computational complexity is equivalent to decoding a single sequence of block length Saej^isCj. However, with the 
sliding window method, decoding complexity can be reduced to C^g (the largest block size) where it is presumed that 
40 the state metric calculations and trace back can be done in parallel so that no additional clock cycles are required to 
conduct the trace back. Note that state metric calculation needs to be done only once for each transition input. Only 
the initial trace back point varies: either state 0 or the maximum state metric state corresponding to either trace back 
from block end or from the final edge of the window; as in Figure 7. 

[0040] Figure 1 a shows an example of hardware implementation of the BTFD. For each input data, 256 state metrics 
45 need to be updated and have to be stored In a RAI\/I with 256 entries. As shown by the butterfly of Figure 3, two previous 
state metrics are necessary to update one state metric. It is obvious that two consecutive state metrics can be updated 
from the same two previous state metrics. As a state metric buffer, one 256-entry RAM is employed. One read and 
one write of previous and new state metric values must be done in a single clock cycle. The minimum and the maximum 
state metrics among the 256 state metrics are calculated and stored to registers. The minimum value obtained from 
50 the previous 256 state metrics stored in the register is subtracted from state metrics read from the RAM. Thus metrics 
at time 1+1 can be expressed as 



metric(i+1 ,2n) = 
max{metric(i,n) -a„,Ji)+Zo[2Yo(n.2n)-1l+Zi[2Yi (n,2n)- 
1]+Z2 [2Y2(n,2n)-1], 
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metric(i,n+1 28) -a„,n(i)+Zo[2Yo(n+1 28.2n)-1 ]+Z^ [ZY^ (n+1 28,2n)- 



1 l+Zg [2Y2(n+1 28,2n)-1 )} 



5 



metric (1+1 , 2n+1 ) = 



max{metric(i.n) -a^,n(l)+ Zo[2Yo(n.2n+1 )-1 ] + [2Y^ (n.2n+1 ).1 ] + 
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Z2[2Y2(n.2n+1)-1]. 



metric(i,n+12B) -a^,n(i)+ Zo[2Yo(n+128, 2n+1)-1] + 



15 



Zi[2Yi(n+128.2n+1)-1] + Z2[2Y2(n+128,2n+1)-1] } 



where af„in(i) is the minimum metric at time i. This subtraction Is beneficial to prevent overflow and to simplify the 
reliability calculation for r (n^rid) • 

[0041 ] The information of path selection for one input data (one transition) requires 256 bits. The sliding window size 
20 and the window overlap size are assumed to be 80 and 64, respectively. 1 6 decoded bits are obtained from a 80-blt 
traced-bacl< sequence. The required size of the trace back buffer is 80x256 bits. With an additional 256-blt storage to 
realize a ring buffer, the size becomes 81x256 bits = 1296x1 6 bits, so use a 1296x1 6 bit RAI\^ for the trace back buffer 
(see Figure 1 a). These additional 256 bits are necessary to provide enough hardware clock margins to complete tracing 
back concurrently. In the worst case, two trace backs will be required from state 0 or the maximum state metric state 
25 at the edge of the window. About 80x2 = 1 60 clock cycles are consumed for the worst case tracing back. The number 
is below the available clock cycles of 256. 

[0042] After state metric calculation for the first 80 received symbols, the trace back buffer is filled up and ready for 
tracing back. After this first trace back, the trace back buffer becomes ready at the completion of every 1 6 transitions 
(256x1 6 state metric calculations). There are two types of trace back operations: (1 ) when the data comes to the edge 

30 of the sliding window (see Figure 9a) and (2) when the transition reaches to one of the block size candidate (see Figure 
9b). For the first case, the trace back starts from the maximum state metric while, for the second case, the trace back 
starts from state 0. Note that 16 decoded bits for the first case is common to all block length candidates. Thus the 
difference in decoded bit sequence is always bounded within the decode sequence in the second case. From this 
reason decoded bits for the first case and the second case are stored in the output RAM separately. Also, the results 

35 of CRC calculations for the first case can be stored and not necessarily be recalculated for each block length candidate. 
[0043] If there are several decoded sequences whose CRC checks Indicate no error, then a sequence with the 
smallest reliability figure is selected. If alt sequences cannot pass the CRC check, then reliability figures of all sequences 
are compared and a sequence with the smallest reliability figure is selected. 

[0044] Figure 1 0 shows the memory map of the output buffer. The decoded bits for the first case are stored between 
40 0 and 27 (area 1 ). The rest of the RAM forms a ping-pong buffer (area 2) so that the most likely (In terms of CRC match 
and smaller reliability figure) decoded bit sequence can be retained In one of the ping-pong buffers. The other buffer 
is used to continue decoding, evaluation of CRC and reliability figure for the rest of block length candidates. Once the 
block size is resolved, how to concatenate decoded sequence from buffer area 1 and area 2 is straightf onward. 



45 



Claims 



1 . A method of determining the length of a received encoded block, comprising: 



55 



50 



(a) receiving a sequence of encoded samples, said sequence Including blocks of a length selected from a set 
of candidate lengths; 

(b) Viterisi decoding said sequence with a sliding window up to a first of said candidate lengths; 

(c) evaluating said decoding of step (b); 

(d) Viterisi decoding said sequence with a sliding window up to a next of said candidate lengths; 

(e) evaluating said decoding of step (d); 

(f) comparing the evaluating of step (e) with prior evaluating; 

(g) repeating steps (d)-(f) to include all of said candidate lengths; and 

(h) determining a block length from said comparing of step (f). 
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2. The method of claim 1 , wherein: 

(a) said evaluating includes checl<ing for errors using parity bits 
5 3. The method of claim 1 , wherein: 

(a) said evaluating includes comparing state metrics at the candidate length. 
4. A method of decoding, comprising: 

10 

(a) providing a set of candidate block lengths; 

(b) providing a window length and a window slide length; 

(c) updating state metrics with subtraction of minimum metric value of previous Interval and starting at state 
0 for a trellis using a sequence of received samples over an interval of length equal to said window length, 

IS said updating including storing corresponding transition bits; 

(d) tracing back from said updated state metrics with subtraction of minimum metric value of previous interval 
using said stored transition bits; 

(e) retaining bits of said tracing back from an end of said tracing back for a length equal to said window slide 
length; 

20 (f) updating said state metrics over a next interval of length equal to said window slide length, said updating 

including storing corresponding transition bits; 

(g) repeating steps (d)-(f) until said updating reaches a first of said set of candidate lengths; 

(h) tracing back from state 0 of said updated state metrics of step (g) using said stored transition bits, and 
combining with said retained bits for a first decoding up to said first of said set of candidate lengths, and 

25 evaluating said first decoding, and continuing updating said state metrics to the end of said next interval of 

length equal to said window length; 

(i) repeating steps (d)-(f) said updating reaches a next of said set of candidate lengths; 

(j) tracing back from state 0 of said updated state metrics of step (i) using said stored transition bits, and 
combining with said retained bits for a next decoding up to said next of said set of candidate lengths, and 
30 evaluating said next decoding, and continuing updating said state metrics to the end of said next interval of 

length equal to said window length; 

(k) repeating steps (i)-(j) ^^^^^ of said set of candidate lengths has been traced back; and 
(I) comparing said evaluatings to select one of said set of candidate lengths. 

35 5. A logarithm approximation method, comprising: 

(a) providing N-bit numbers c and d where N a positive Integer and d is less than or equal to c; 

(b) extracting M-bit number m from said c where M is a positive integer less than N , said m aligned to the most 
significant non-zero bit of c; 

40 (c) extracting M-bit number n from said d, said n aligned to said m; 

(d) piecewise linearly approximating a logarithm of said m; 

(e) adjusting the approximation of step (d) by a table entry; 

(f) piecewise linearly approximating a logarithm of said n; 

(g) adjusting the approximation of step (f) by a table entry; and 

(h) approximating a logarithm of the quotient d/c by the difference of results of step (g) and step (e). 

6. The method of claim 5, wherein: 

(a) iVl equals 5; and 
so (b) said logarithm has base such that 1 0logi^332 = 16 

7. A method of logarithm approximation, comprising: 

(a) providing N-bit positive integer c = x.,X2 . . . with N a positive integer and x^ the most significant bit and 
ss Xm the least significant bit; 

(b) providing N-bit positive integer d = yiy2.. with y^ the most significant bit and y^ the least significant bit 
and with d less than or equal to c; 

(c) extracting l\/)-bit number m from said c, said m with most significant bit aligned to Xj of said c where x^. 
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X2, x,^.^ are all 0 and is 1 where k is a positive integer with 1 < k ^ N; 

(d) extracting M-bit nunnber n from said d, said n with most significant bit aligned to yj^ with said kfrom step (c); 

(e) when said m Is in the range 2> to 2'+^ -1 for a positive integer I, define a number t = m2W-i + (M-i)2M"^ + 
table [m] where table[ ] Is a lookup table; 

5 (f) when said n Is In the range 2l to 2i+^ -1 for a positive Integer J, define a number u = n2M-i-j + (j-1)2W-i + 

table[n]; and 

(g) taking t-u as an approximation of log(d/c). 
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Figure 1a 
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Figure 1 b 
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Figure 2 
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Figure 3 
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Figure 4 
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Figure .S' Basic processing flow of blind transport fonnat detection 
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Figure 9a 
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Figure 10 
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